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^ : Abstract 

\^ • We derive the hard thermal loop action for soft electromagnetic fields in the 

finite temperature world-line formulation at imaginary time, by first integrat- 
^ , ing out the hard fermion modes from the microscopic QED action. Further, 

using the finite T world-line method, we calculate all static higher order 
' terms in the soft electromagnetic field. At high T, the leading non-linear 

■ terms are independent of the temperature and, except for a term quartic in 

the time component of the vector potential, they cancel exactly against the 
, vacuum contribution. The remaining T-dependent non-linear terms become 

r"! \ more strongly suppressed by the temperature as the number of soft fields in- 

creases, thus making the expansion reliable. Applications of this method to 
other theories and problems at the soft scale are also briefiy discussed. 



.!-H . I. INTRODUCTION 

X 

d \ An important feature of interacting field theories at finite temperature is the generation 

of thermally induced mass scales, proportional to the coupling constant g When an 

ensemble of interacting particles at temperature T is disturbed, the system responds by 
creating excitations at the energy scale gT <^ T. Here we assume g <^ 1. Consequently, 
there appears another mass scale due to interactions, the soft scale gT, in addition to the 
hard scale T that characterizes the average energy of the constituent particlesQ Since {gT)~^ 
is much larger than the average particle separation T~^, it is clear that such an excitation 
is a collective phenomenon of the system. 

A natural question then is how to describe the physics at the soft scale in terms of 
the effective degrees of freedom. This problem has been studied intensively during the last 
decade, starting with the identification of the relevant Feynman loop diagrams that has 



^This statement is true for QED. In general, there may be additional scales that arise due to 
interactions. In QCD, one also has scale g^T due to screening of magnetic interactions. 
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to be taken into account in amplitudes involving soft external fields Named Hard 

Thermal Loops (HTL), the diagrammatic approach then led to considerations of the low 
energy effective action that captures the physics at the soft scale. 

Such an effective action was first considered in , and a very simple expression was later 
obtained in [^. For example, the behavior of the soft electromagnetic field is described by 
the following gauge invariant effective Lagrangian, 



^soft 




(1) 



where F^^ is the field strength tensor, = (z, k) with k the unit three vector, the partial 
derivative and mj^ = the thermally induced Debye mass. The HTL action was later 

on derived in a variety of approaches, for example, from transport theory by truncating the 
Schwinger-Dyson equations P as well as by purely classical kinetic treatments 0. 

In this paper, we demonstrate that the HTL effective theory can be obtained rather 
straightforwardly from the full microscopic action by simply imposing a momentum cut-off 
and integrating out the hard modes. Such a procedure is close to the renormalization group 
picture of Wilson (at the one- loop level), although we will neglect all cut-off dependent 
terms that are subleading compared to T in the resulting determinant. In this respect, our 



approach is similar to the derivation in |]T0[. However, we use a completely different method 
to evaluate the determinant, namely the world-line formulation, wherein the effective action 
is written quantum-mechanical path integral. 

At T = 0, the world-line method has been developed and successfully applied to a variety 
of problems for a long time (see []TT| and references therein). In [|T^]-|T^ the formulation was 



extended to finite temperatures, and the finite T method was subsequently used in several 
studies of the effective action [|TR|TB[. One of our main objectives here is to demonstrate 



that this formalism has further feasible and attractive applications at finite temperature. 
As we show with the example of the HTL effective action in QED, the finite T imaginary 
time world-line method can be useful even for calculations of non-static quantities. 

In addition to the HTL effective action, we also discuss the higher order, non-linear terms. 
Indeed, one of the true advantages of the world-line method, compared to diagrammatic 
calculations, is that it comprises the sum of all 1-loop diagrams to a given order in the 
external field. Using the relative simplicity of the world-line method, we demonstrate how 
the static contribution to the effective action, for an arbitrary number of soft external 
photons, can be obtained. Although these terms beyond quadratic order are well known 
to not be of the HTL type, the non-linear interactions are needed for studies of light-light 
scattering and photon splitting (see fl^ and references therein). Such non-linear processes 



have by now been observed in laboratory experiments |jT8|| , and they may also be important to 



explain astrophysical effects such as the observed 7-ray spectrum from very massive objects 
||19|| . Even though the kinematical region considered here may not be encountered around, 
for example, neutron stars, it is nevertheless interesting to explore the form of the effective 
action under different circumstances and obtain a more general and complete picture. In 
that sense, our results reported in this paper adds rather nicely to the earlier investigations 

[Ml- 



m 



At this point we should mention that we have not been able to generalize the non-static 
case beyond the hard thermal loop action. Including the dynamical information is clearly 
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of outmost importance. Hopefully the results reported here can serve as a first step in that 
direction. 

Finally, it should also be noted that the world-line formalism provides an intuitive con- 
nection between classical kinetic theory and the underlying quantum field theory. At 



finite T, this was noted in ||2^ and later explored in |^. Remarkably, a world-line ap- 
proach, analogous to that used to obtain classical kinetic theory, can be applied to derive 
an effective action for small-x physics ||26|| . 



The paper is organized as follows. In the next section, we recall how to go from determi- 
nants to quantum-mechanical path integrals and write down the one-loop effective action 
for soft gauge fields in finite temperature QED. In section III, we explicitly derive the HTL 
effective action from the one loop effective action. In section IV, we go beyond the HTL ap- 
proximation and obtain a general expression for the effective action to all orders in the static 
component of the gauge field. We end in section V with our conclusions and a discussion of 
further applications. Details of the computations are discussed in the appendices. 

II. FROM DETERMINANTS TO QUANTUM MECHANICAL PATH 

INTEGRALS 

The usefulness of rewriting one-loop effective actions in the language of first-quantized, 
quantum-mechanical path integrals was first emphasized by Strassler He showed that 
the method was equivalent to the Bern-Kosower rules for one-loop n-point scattering 
amplitudes in gauge theories, the latter having been derived originally within the framework 
of string theory. In this section, we will begin by briefly reviewing the work of Strassler. For 



a comprehensive review of later developments, we refer the reader to |TT|]. Following p2| , p!3 
we then proceed to the finite temperature imaginary time case. 

A. Scalar QED at T = 

To be specific, consider as an example scalar QED in four Euclidean dimensions. The 
effective action Feg for an electromagnetic background field is. 



eff 



log 



(2) 



where = + with the electromagnetic vector field and e the electric charge. 

Performing the Gaussian integration, and using the Schwinger proper time method, we 
obtain 

Feff = logDet \-D^] = - [d^x r -(a;|e-*(-^')|a;) . (3) 

L -I J Jot 

We can interpret the matrix element in Eq. as a transition amplitude, and express it in 
the form of a quantum mechanical path integral 



Teff = - y / d^x U Vx e-^->- = - j Af J^^^,^^ Vx e"^" '- , (4) 
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with the world-hne action S'w.l. given by the expression 



•S'w.l. = / dr 
Jo 



'^^ — iex^r) ■ A{x{t)) 



(5) 



Here = dx^/dr. The normahzation constant A/" in the effective action ensures that the 
free (e = 0) theory satisfies the condition 



Af 



2.V(i)= 
x>(0)= 



Vx e 



1 



(6) 



Physically, Eq. (^) can be interpreted as the path integral for the first-quantized theory of 
a particle interacting with a classical background field. The endpoint positions of the particle 
are constrained by a;(0) = x{t) = x, with an integration over all possible locations x and 
intervals t = At. It should be noted though that the variable r does not represent any real 
time but is merely a parametrization of the trajectory. As such, a re-parametrization r — »• 
f(r) does not affect the classical action if the associated square-root of the one-dimensional 
metric, the einbein e, transforms in the proper way. However, the complete path integral is 
not invariant unless one integrates over the space of all metrics as well |]3^. In the case at 
hand, the einbein has been fixed, without loss of generality, to e = 2. 

The above method can easily be generalized to include spinning particles as well. This is 
achieved by representing the trace over spinor indices as a path integral over anticommuting 



This procedure 



Grassmann fields in a supersymmetric generalization of Eq. (^) 
applies in general to internal symmetries, and as was shown in [^, even additional internal 
indices such as color can be represented by a true scalar world-line action. 



B. QED at Finite Temperature 

We will now derive the corresponding form of the effective action at finite T for soft 
external photons in QED, where the soft photons have a momenta p ~ eT ^ T. We will 
further assume that the soft momenta are much larger than the electron mass, p ^ m, so 
that the electrons can be considered massless. The microscopic, finite temperature imaginary 
time action is 



'5'micro — J dx^ I d X 



(7) 



where X4 G [0,/5 = 1/T] and F^^, = d^A^ — d^A^, with the tildes for later notational sim- 
plicity. As is well known in the imaginary time formalism, the gauge fields are periodic in 
the compact direction, y4^(x, X4 + (3) = y4^(x, 0:4), whereas the fermions are anti-periodic, 
X4 + jS) = -^/'(f , X4). 

Since both the fermion and gauge fields contain all possible momenta, we pick a scale A, 
such that eT -C A -C T, and then split the fields into two parts; one part contains the soft 
momenta and the other, to be integrated out, contains the hard momenta. By performing 
this cut-off procedure, we introduce an explicit A-dependence into the effective action and 
also break the gauge invariance. However, the terms in the resulting effective action that 
depend on the scale A are subleading, since A/T ^ 1. They can therefore be ignored here 
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since we will only be interested in the leading dependence on T. The remaining terms, being 
independent of A, thus have to be gauge invariant. 

To derive the one-loop effective action for a soft electromagnetic field one can neglect 
the hard photons. We therefore take A^ = A^. Simultaneously, we write %p = ip + , where 
\E' contains the hard momenta. We then use the background field method to expand 
the action in Eq. (|^ to quadratic order in Integrating out the hard field, we obtain the 
effective action for the soft photon modes 

S = Jd^x Qf^^ + #7/.(5m + ^eA^)^^ - log Det [tip] = S^^ft + ^soft , (8) 

where Ip = 7^(5^ + ieA^) and ^soft contains the determinant. The determinant can be 
re-written as 



Det [tip] = V(Det [tip])' = ^Det [-DH + tea^.F^,] = ^Det [-D^] {l-tea^,F^jD^) , (9) 

where cr^^ = [7^,7^]/4. 

Now, acting on the hard field gives a contribution ~ T, whereas the magnitude of 
the soft field A^ is ~ T. By power counting arguments, it follows that the term eF^y/ D"^ is 
sub-leading compared to the 1 in Eq. (|]) when it comes to the HTL effective action. Thus 
the dominant correction from integrating out the hard fermion field becomes (after taking 
a trivial trace over spinor indices) 

S'soft = -2Trlog[-Z}2] _ (IQ) 

Apart from an overall numerical factor, the important difference between Eq. ([T0| ) and the 
corresponding equation at T = 0, Eq. (^), stems from the fact that the finite T determinant 
has to respect the appropriate periodicity conditions in the compact X4-direction. In this 
case, it therefore has to be evaluated on the space of anti-periodic functions, due to the anti- 
periodic boundary condition of the fermion field. As was shown in ref. [1T3, the resulting 
quantum-mechanical path integral is 

Ssoft = 2 E i-^r / -AT/ Vxe-'-- , (11) 

n=-oo 1^ •'x^{0)+n/35^4 

where S'„.i. is the world-line action in Eq. (|). 

In the imaginary time formalism of world-lines, a finite temperature corresponds to the 
fact that the path can wind an arbitrary number of times around the cylinder x before 
returning to the starting point. For n = there is no winding, so that particular term is 
insensitive to the compact form of 0:4 and therefore gives the T = contribution. 



III. HARD THERMAL LOOPS FROM WORLD-LINES 

In this section, we will derive the HTL effective action for soft photons from Eq. ([TT]). 
Many of the details of the computation are given in Appendix A. We begin by noting that 
Eq. (|l^) is the expectation value for the Wilson line, and can be expanded in powers of A^, 
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5'soft = (e*''-^o "^y = (5'soft(o) + 'S'soft(i) + 5'soft(2) + • • •) , (12) 



where the subscript [i) indicates the number of external fields. The individual contributions 
are all gauge invariant under the periodic gauge transformations A^ — {1/ e)d^uj, where 

u!{x, X4 + n(3) = uj{x, X4). 

To lowest order in the coupling constant, the path integral describes a free theory, and 
the corresponding contribution to ^soft is easily obtained. Performing the substitution 



(r) + np6^iT/t + , (13) 



with Zfj^ a r-independent constant and u^{t) = u^{0) = 0, we use Eq. (|) to get 

^E(-l)"f f^-"'"'^'. (14) 



V4 p dt 
H^^ 8^ Jo 



with V4 the four- volume. The first term is ultraviolet divergent but independent of the 
temperature, whereas the second term gives 



5'soft (0) 



T>o ~ 180 



(15) 



which is the expected result for free fermions when the infrared cutoff A is neglected. 

The next term 5'soft (1) in Eq. (|12|) is linear in the soft external field and vanishes, due the 
sum over winding modes being odd and /q driu^^Ti) = u^{t) — u^{0) = 0. This result is of 
course just a manifestation of Furry's theorem. 

Now consider the term quadratic in the external field, S'soft(2)- For simplicity, we will 
work in a gauge where ^4 is time independent, namely ^4^4 = 0. The spatial part Ai 
[i = 1, 2, 3) of the vector field remains nonstatic and transforms under the residual static 
gauge transformations as Ai ^ Ai — (l/e)Via;. There are three kinds of terms to consider: 
(A4-A4), (Ai-Aj) and (^4-^4). These will be discussed below, with some of the details of 
the derivation given in Appendix A. 



A. The ^4-^4 term 



Expanding out the effective action to quadratic order, the most general form of the 
effective action is given by Eq. (A. 2). For the term containing two powers of A4, Eq. (A. 2) 
can be simplified to read 

SsoM2){A,-A,) = -2e' f f dz,A^^ A^ £ (-1)" H ^e-'^^'^Ar f dr, H dr, x 

J (ZTCy Jo Jo t Jo Jo 



j Vuiu^{Ti) + '^\ L4(r2) + ^)exp - j'^ dr [u" /A - ip ■ u{5{t - T2) - 5{t - n))] 



n=— 00 
t 



(16) 
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where A^P = A^{p} and Af^ = A^i-p). 

Following Strassler |5S|, we now exponentiate the vector potential under the condition 
that only terms linear in each ^ contribute: 



ujTi) + 



A 



,(«4(rO + ^)A« 



linear in each A 



(i) 



(17) 



As discussed in appendix A (see Eqs. (A.4)-(A.10)), this exponentiation allows one to re- 
write the exponential in Eq. (|T6|) in terms of one dimensional Green's functions of a free 
particle constrained to a circle of circumference t. Subsequently, the path integral can be 
made Gaussian after a change of variables. Since the Green's function and its derivatives 
are known, one finds, by defining 7 = p'^x{l — x) and using Eqs. (A. 13) and (A. 14), that 
the finite temperature part (n 7^ 0) of the effective action (|16D becomes. 



soft (2) 



(A4-A4 



T>0 



00 

n=l 



(18) 



use 



The leading contribution in Eq. (18) is obtained when (3^/^ 0. In this limit, we can 

2 1 



lim K2{n(3^) = 



lim Ki{n(3y/y) = „ ^ , 



(19) 



since the sum over the winding modes remains convergent. Employing the identity, Jq dx{l — 
x)[l — 6{x)] = 0, with the usual convention [p!2| , p8| for the sign function, e(0) = 0, we obtain 
finally the result 



5'soft{2)(^4-^4y 



T>0 



6 



(27r)3 



n=l 



Ai{p)Ai{-p) dz4 = mjj J d^z[A4{z)Y 



(20) 



where mj^ = is the Debye mass. In the static limit, this result is well known to be 

the only leading contribution and gives rise to a screening of the potential between 

two static test charges. 



B. The Ak-Aj term 

The derivation of this part of the HTL effective action is not as straightforward as the A4- 
A4 term because the spatial components of the vector field also depend on x^, Ai = Ai{x, X4). 
Using Eqs. (A. 12) and (A. 15) in Appendix A, the effective Lagrangian for this term becomes 

tr^ , I Lib „2^2/A-t I , . -ip4,nl3x^-tP'^x{l-x) 



L.oM2Mk-A,) = E dx{l-x)e- 

x|pfcp,4^)Af [l-2a;]2 + ^l«.i(2)[l-5(x)]} , (21) 
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where the shorthand notation is A^^^ = A{P) and A^"^^ = A{—P), with a Euchdean 
four- vector, = pl+ p"^ and = uJi = 2'kIT the external Matsubara frequency. 

We now spht the effective Lagrangian (|2l| ) into a piece L^^-* for the term including the 
delta function, and another L^"^^ containing the rest of the expression. The finite T part of 
L^^^ can be evaluated directly: 



-(1) 

^soft (2) 



— e 



1(1). 1(2) g (_l)n ^g-n2/3V4i 

n=-oo ^ ■'^ 

1(1). 1(2) _ 



-X e 



-ip4n/3x— tP^x(l— x) 



12 



T>0 

(22) 



This term looks like a magnetic mass-term, but since such a term is forbidden in QED ||3J] it 
has to be cancelled by the contributions from L*-^-*. In addition, it also violates the residual 
gauge invariance. 

As in [13 1, the remaining terms in 

L(2) 

can be simplified further by first using the 
substitutions n —n and y = 1 — x. If f{x,n) denotes the part of L^"^^ that only 
has to be summed over n and integrated over x to produce the Lagrangian, we have 
J2nIodx{l — x)f{x,n) = J2n lo dyy f {y , n) , by using the fact that p4 is a discrete Mat- 
subara frequency. As in the purely static case, we can then write the result in terms of 
Bessel functions of the second kind, to find the finite T contribution 



-(2) 

^soft (2) 



X 



r dx{l - 2a;) £(-l)"iro (n/3y'P%(l 

•^0 n=l ^ 



X] 



PkPjA^f^'^Af\l - 2x) cos[npiPx] + 2p2l(i).l(2) £ da cos[npiPa]^ . (23) 

To obtain the second term in the above equation, we performed an integration by parts in 
X, and then defined the integral over the auxiliary variable a to re-write a sine-function as 
a cosine. 



We now use the following identity for the infinite sum over winding modes |35 

1 



^ (-1)" cos[nK\Ko{n(, 



n=l 



^ m=0 



2 ^ 

^ m=0 



1 



^(l)^ + [{2m + l)n-K] 
1 



+ 



+ [(2m + l)7r + K] 



1 OO 

-E 

^ m=0 



[m+ 1) 



(24) 



where 7^; is Euler's constant. In fact, this identity converts the sum over winding modes into 
a sum over fermion Matsubara frequences. Keeping only the leading temperature dependence 
of Eq. (pW), we show in Eqs. (A.16)-(A.19) that the finite T Lagrangian (53) becomes 



-(2) 

^soft (2) 
/ „2 



T>0 



2^2 A^ 



^2 



:T E 

3 m=— 00 



00 rl 
2 



dx 



x(l — x) 



log P^x(l -x) + e'^ + (ki - xpif 



(25) 
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where = (2m + l)7rT plays the role of the fermion Matsubara frequency. 

Now, to extract the HTL effective action from Eq. (pSD we note that the leading con- 
tribution is equivalent to a hard thermal loop integral, multiplied by A^j}^A^^\ This can be 
seen as follows. First take the derivatives with respect to pk and pj to get 



-(2) 

^soft (2) 



T>0 



.(1) .(2) 



T E 



m=— oo 



dee^ 



dx 



kj 



P'^x{l - x) + 6"^ + {ki - xp^y 



2p.jPkx{l - x) 



[P^x{l 



+ (^4 - xpif] 



Partially integrating over 9 in the first term in Eq. 
independent UV-divergent constant, we are left with 



(26) 



q), and neglecting a temperature 



-(2) 

^soft (2) 



T>0 



if! 

37r2 



m=— oo 



dee' 



dx 



e'^Skj - 3pjPkx{i 



XI 



[P^x{l - x) + e^ + {k4 - xp^) 



212 



(27) 



Since the ^-dependent integrand only contains multiples of 6'^, we can formally interpret 
e as the magnitude of a three- vector, so that e"^ |6'p and dee"^ d^e/ (47r). By adding the 
terms A^f}^ Af\x{ekPj + ejpk) — ekPj], which are all odd in e, and changing the integration 
variable e = k — xp, we obtain the leading contribution. 



-(2) 

^soft (2) 



T>0 



oo „ 



d^k 



kkkj 



(27r)3 K^{K - Py 



(28) 



where = k\ + k"^. Once the sum has been performed, one can continue p^ to arbitrary 
Euclidean values, so that it becomes the analytical continuation of the energy in Minkowski 
space-time. When both p4 <C T and \p\ <C T, as is the case when the fields are soft, the 
integral in Eq. (pSD is well known p|. The leading contribution from Eq. (|25|) is 



-(2) 

^soft (2) 



T>0 



ipi 



K -P 



k ■ A{-P) . (29) 



Thus, as anticipated, the first term in (^) cancels the result in Eq. (|2^). The total 
effective Lagrangian for the spatial part of the electromagnetic field is then 



^soft (2) 



{Ak-Aj 



T>0 



/ — k ■ A(P) 



jp^ 

k-p 



k-A{-P) . 



(30) 



This term is gauge invariant under the residual static gauge transformations allowed by the 
condition 94^44 = 0. 



C. The A4-Ak term 

Finally, we have to consider the part of the effective action that contains one spatial 
component Ak together with one power of A4. Since the y44-field is static, we have after a 
Fourier transformation (see Eq. (A. 2)) 
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(2.V n 

n(3 



(31) 



However, the n/3/t-term vanishes since the sum is odd in n, whereas the 'it4-term gives, after 
an exponentiation in the same way as before. 



'S'soft(2)(^4-^fe) 



71^ J (27r)3^' ' Jo 



t /-n 



JO 



linear in 
Afc and 



(32) 



But this part of the effective action has no dependence on ^44. Therefore, the exponent 
cannot be expanded to give any term hnear in ^4, and the Lagrangian vanishes to all orders 
in p/T. 



D. Summing the contributions 

The complete effective action to quadratic order in the soft photons then becomes, by 
adding the contributions from Eqs. (pO[) and (|30| ) 

5soft = mlj d'x (^[A,{x)]' + J^k- A{x) j-^^j k ■ A{x)^ . (33) 

We have thus derived the effective action by first integrating out the hard modes from the 
microscopic action, and then evaluating the determinant form of the one-loop effective action 
in Eq. (H) as a quantum-mechanical path integral. The general HTL action in Eq. (P of 
course reduces to the above result when the gauge condition 94^44 = is imposed. 



IV. BEYOND HARD THERMAL LOOPS 

In this section, we turn to the higher order terms in the effective action. Although these 
terms do not have HTL's, it is nevertheless interesting to study the form of the effective action 
in the high T phase, and obtain a systematic expansion in powers of p/T. As mentioned 
already in the introduction, the non-linear effective action is a useful tool in astrophysical 
applications pO|]. Moreover, one could use the resulting non-linear expression to contract 



and thermalize any pair of external legs, and by this procedure obtain higher corrections to, 
for example, the pressure. 

The power of the world-line approach becomes manifest as one goes to higher orders. 
For instance, to obtain the dominant term at high temperature, one should note that indi- 
vidual loops in the diagrammatic approach may contain superficially leading terms that only 
cancel when all contributing graphs are summed. In marked contrast, within the world-line 
formalism, all diagrams are accounted for in the effective action. 
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In [0,0] a diagrammatic approach was used to argue that all terms with > 4 external 
photons have a finite limit when T ^ oo. In the world-line formulation, we will verify 
explicitly that this is true for the static terms. In principle^], we can also obtain all corrections 
in a power expansion of p/T. For iV = 4, parts of our results overlap with the world-line 
calculation in 



16|. 



Proceeding in the same way as before, by generalizing straightforwardly Eq. (A. 12) of 
Appendix A to the A^'th order, we have in the static case 



-te 



soft (N) 



dxA 



d^Pi d^pN 
(27r)3'"(27r)3 



(53(j9i + ...+PN) E(-l)" / dtt^'-^e-''"'''/''' X 



exp 



^ i=l 



dxy ■ 



XM-l 



dxjy exp 



N 



{tpi- Pj{Xi- Xj){l - 



^.))} 



i,j=l 
i<j 



Kj 



6{xi - Xj)) 



(34) 



where A^^^ = A(pi), and the only term in the field-dependent exponential to be kept is 
the one linear in each A^^\ Since the effective action refiects the symmetries of the original 
microscopic action, and in particular the invariance under charge conjugation, all odd powers 
of the vector field vanish and the only contribution comes from A^ even. This generalizes 
the explicit calculation for a single external field discussed in section III. 

We now expand the exponential exp(n/3 X^ili ^i*"*/^) powers of A^. For odd powers, 
the sum over n is odd and vanishes. Hence we only have to consider even powers. The 
effective action could still contain an arbitrary even number k [0 < k < N) of powers of A4, 
though. Following the general procedure outlined in appendix A, we then have 



soft [N) 



57r^ Jo 



dxA 



(2^---(2^^^^^ + ---+^"Vo^"^--7o 



N 



dx^Y^^N^ (35) 



k=0 
fceven 



where the function consists of three different terms. 



Ak) , Ak) , Ak) 



E 



Ak,n) Ak,n) Ak,n) 



(36) 



The three terms in Eq. (^^ 
tial in the last line of Eq. 
exp{nj3J2iLi ^4^^ /t) is expanded to an arbitrary even order k. When only the first p- 



correspond to the different ways of expanding the exponen- 
to get A^ external fields in total, given that the term 



dependent term of the last exponential in Eq. 
powers of the vector field, one finds. 



contributes with the remaining (A^ — k) 



^With the operator D^, this is true for scalar QED whence the boundary conditions in X4 is 
changed. In QED, one would also have to take into account the -F^,y-term to get the correct factor. 
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(37) 



and if the {N — k) powers only come from the p-independent term of the exponential, the 
result is 



4^-]= r ^^-^•^/^(r^/5)^^+'='/^-^i^2-.(.-^)/2 {nm Mix}) 
Finally, there is also the possibility to expand both terms, giving 

(N-k)/2-l 



(38) 
(39) 



/Jfg;=e[(iV-4)-A:] P 

1=1 



n(2+'=+'-^)(n/3)^-'-^ir2+fc+._^ {n(3p) ff{{x}) x 



V 



(40) 



where the ^-function ensures that the term vanishes when k > {N — A). In all the equations 
above we used 



N 



P"^ = {Pi-Pj^ 



Xj))} 



(41) 



i<j 



and the functions fi{{x}) comprise the dependence on the dimensionless variables Xi, to- 
gether with some irrelevant combinatorical factors. V denotes both all permutations of 
choosing k factors of ^4 and the other (A^ — k) factors out of the total, as well as all pos- 
sible ways of contracting the different vector fields with themselves and/or the momentum. 

To study the high temperature behavior of the effective action, we isolate the depen- 
dence on T and p in From Eqs. (p7D-(^) we find that all three terms give different 
contributions. As shown in appendix B, the p and T dependence of -^^(1) becomes 



1%) ^ fi{{x}){A,,. . .A,J{A^^^. . -A, 



.P. 



Af-4 



ioo+e 



2m 



zT{z){l - 2'-^)az) {-z-l)...{-z-N + 3) {[3pr x 



JkN 



Af-4 



2T[{N-4 + z)/2] T[{N + z-2)/2] 



+ h{N-2)PiN^^PiNp 



T[iN + z-2)/2] 



T[{N + z-3)/2] T[{N + z-l)/2]_ 

W + z)/2] 



j=0 
j even 



r[(j + ^ + i)/2] 



E 

m=l 



(1). 



r[(iv + z-i)/2] 

r[(j + Z + 2m)/ 2] 



r[(j + ^ + 2m + l)/2] 



(42) 



where C,{z) is the Riemann zeta- function, a^-* contains some irrelevant combinatorical factors 
and the sum over m runs up to the largest integer < [N — j — 2)/2. Similarly, one can show 
that the expression for the second term /^(2)' valid for all if > 6 and for k > 2 when 
= 4, can be written as 
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Af-4 



—zT{z){l - 2'-^)az) i-z-l)--- i-z -{N + k)/2 + 3) x 

— ioo+e ZTTl 



>kN 



N-4 



2T[{N-A + z)/2] r[(iV + ^-2)/2] 
T[{N + z-3)/2] r[(iV + z - l)/2] 



r[(iV + z-4)/2] , 
+ ^'=(^-2)r[(iV + z-3)/2] + 



J even 



r[(j + z)/2] r[(j + z + 2m)/2] 



T[{j+z + l)/2] 



m=l 



T[{j + z + 2m+l)/2] 



(43) 



When = 4 and /c = it is easier to sum the Bessel function K(j(n(3p) in Eq. ( p9| ) directly. 
The leading high temperature term contains a logarithmic dependence on T, but when 
integrated over Xi this contribution vanishes. This result, with four external fields, agrees 
with the result in |T^. Finally, the last term /^jg) becomes, for a given k < N — 4, 

(Af-fc)/2-l 



1=1 

2{N-l)-i-k 



27ri 



^r(^)(l - 2i-^)C(^) (-Z - 1) ■ ■ ■ {-z - (iV - /) + 3) (/5p)-^ X 



V[{k + z)l2] , (^-'^^-2)/2_^(3) r[(A; + ^ + 2m)/2] 



r[(A; + 2+l)/2] 



E 

m=l 



r[(A; + z + 2m+ l)/2] 



(44) 



Even though the Eqs. (^21)- (^) look complicated, the analytical structure of the integrals 
allow some simple but powerful conclusions. This is due to the fact that the temperature 



dependence is only contained in the term {f3p)~^, and in each of the terms in Eqs. (^2]) - 
(0) the contour can be closed in the left side of the complex z-plane and evaluated by the 
residue theorem. The analytical structure of singularities in the complex 2;-plane therefore 
determines the temperature dependence of 5'soft(Ar)- 

The common function g{z) = zr{z){l — 2^^^)({z) has simple poles aX z = ~{2n + 1), 
with n any non-negative integer. Further, the other F-functions that enter the expressions 
always appear in specific ratios, so that the zeros of the denominators in these ratios exactly 
cancel the poles of g{z) when the F-function in the numerator starts to have poles along the 
negative real axis. This holds irrespective of the additional factors [{—z — l){—z — 2) ■ ■ ■] in 
the numerators that actually cancel several of the poles. Consequently, there are only simple 
poles to evaluate. If there were any higher order poles, the factor {j3p)^^ would contribute 
with a logarithmic factor, but since this is not the case, the effective action contains only 
powers of p/T. 

Although the above equations give the complete, general formula, a simpler result can be 
used for the terms independent of the external momenta. This limit corresponds to taking 
the vector field to be constant, and we then have for the effective action, 



S 



2E(-i: 



°o dt 



Af 



2T^Va 



E 



^ienf3A4 



(45) 
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giving for the finite T effective potential l^(A4) 

-1) 



/|T-l4 oo 

E 



IT 



n=l 



■ cos(e?T,/5y44) 



(46) 



Thus the potential is periodic in A^, with period 27rT/e, as was first noted in Restricting 
A4 to the interval [— vrT/e, vrT/e] and writing A4 = nTq/e, — 1 > g < 1, we have 



V{A, 



1 

45 



1 

24 



(47) 



where Qm = {q + l)mod2 (see for a discussion on this effective potential). Expanding the 
potential gives, 



180 



V{A4 



m]j{Ai 



12-k' 



:{A, 



(48) 



This result agrees with the dimensionally reduced theory in ||2^0, once the variables are 
properly rescaled to the canonical mass dimensions in 4d. The result for ^(^4) can of 
course be obtained from Eqs. (^)-(|^) as well, and we have checked that this is indeed the 
case. More importantly though, they also give all the additional derivative terms. 

Apart from the quartic term in Eq. (^8]), there is an infinite sequence, with different A^, 
of leading terms that are independent of the temperature when p/T 0. From Eqs. (^2]) - 
(0) the temperature independent terms must come from the poles at z = 0, corresponding 
to the = terms. Since the = terms are the only ones that also appear at T = 0, 
it is interesting to compare the two contributions. We find from computing the T = 
contribution of Eq. (pS]), arising from the n = term in Eq. (|36|), that all T-independent 
terms in Eqs. (^)-(|4^) are completely cancelled. 

It is straightforward to also include the F^^ term from Eq. (P) into this calculation. In 
the path integral formalism this term will augment the world-line form of the effective 
action by a factor proportional to 



Tr exp 



TV 



(49) 



where k = 1,2,3 and the trace is over spinor indices. Now, this term occurs with the same 
powers of the variables A^, p, t and n as the second term in the last exponential of Eq. (p4D, 
namely the term 



exp 



TV 

{i{p, ■ - pj ■ 1«)(1 - 2{xi - X,))} 



(50) 



The structural dependence on these variables is the crucial ingredient to prove that a can- 
cellation indeed occurs. Since we know that all the T-independent terms in Eq. (R3) are 



For QCD, results for the one loop static potential were first derived by Nadkarni |38| 
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cancelled against the T = effective action, the identical structures of Eq. (^91) and Eq. (|50D 
ensure that this holds true also for the F^^ part. Consequently, apart from the quartic term 
in Eq. (|^), all other higher order terms in the effective action for soft photons in QED are 
suppressed by factors of T. 

As can be seen from Eqs. (^21)- (PD, the polynomial in [{—z — l){—z — 2) ■ ■ ■] cancels all of 
the poles up to order A^, with the exception of the pole at z = 0. Naively, one could imagine 
that the high temperature expansion would yield a series like (1 + + j3^p^ + ...), where 
the first term corresponds to the leading temperature independent part and denotes some 
generic combination of the external momenta. Instead, the leading, temperature dependent, 
contribution for external fields is suppressed by where a is of the order of A^. As 
a result, the terms in the effective action tend to be more strongly suppressed by powers 
of the inverse temperature as the number of external fields, A^, increases. Since the T- 
independent terms by dimensional arguments have the schematic form (9^^~^''(eA)^, they 
are all of the same order for the soft fields, where 9 ~ eT and A ~ T: d'^^~^\eA)^ ~ 
(eT)^^~^)(eT)^ ~ (eT)"^. However, the surviving terms that are suppressed by powers of 
the inverse temperature become smaller in magnitude as the number of the soft, external 
fields increases. The expansion in powers of the soft field is therefore reliable, with the higher 
orders becoming more and more negligible. 

In a cancellation of all non-linear terms was verified in the long wave-length limit, 
IpI ^ 0, 7^ 0- For slowly varying fields, a similar result was found by expanding in powers 
of m/T [|20|, m being the electron mass. We arrive at the same conclusion in the static 
ultrarelativistic limit, apart from the appearance of an (yl4)^-term, by a direct calculation of 
all A^ > 4 terms. This calculation complements the earlier ones and strongly suggests that 
the T independent non-linear terms beyond quartic order always cancel out, regardless of 
the momenta of the external fields. Whether or not the remaining terms are more strongly 
suppressed by the temperature, as in the static case, is not clear at the moment, although 
we expect this to be the case due to gauge invariance. 

V. CONCLUSIONS 

In this paper, we derived the hard thermal loop effective action for a soft electromag- 
netic field directly from the microscopic, one loop effective action in QED. This result was 
obtained by first integrating out the hard fermions from the original microscopic action, and 
subsequently computing the resulting determinant in the quantum mechanical world-line 
formalism. In addition, we have shown how all leading higher order static terms can be 
obtained with the same method. When T — >■ oo, all the remaining terms, with the sole 
exception of the quartic term in the A4-potential, are cancelled by the T = terms in the 
effective action. The remaining terms are generally suppressed by a factor ~ Z^*^, where a is 
of the order of the number of external fields. We would now like to comment on our results 
and speculate about some further applications wherein the methods described in this paper 
may prove useful. 

First of all, we have neglected all dependence on the momentum cut-off, A. This is 
reasonable as long as we are only interested in the leading behavior. In general, the effective 
action should of course also depend on the cut-off that separates the hard modes from the soft 
ones. A simple momentum cut-off does not respect gauge invariance but, in principle, one 
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could derive the effective action in a gauge wliere tfiere are only physical degrees of freedom 
present. In this case, the momentum cut-off does not pose any problems and one would have 
a true Wilson effective action S{T, A) for the soft electromagnetic fields in the QED plasma, 
albeit valid only in that physical gauge. Such a formulation of the HTL effective action could 
be advantageous both for non-static quantities as well as for computations of the pressure 
and entropy in an HTL resummation scheme [^,0. Furthermore, for a quantity like the 



pressure, even the static higher order terms can be used to obtain sub-leading contributions 



41,42 



An interesting point to consider further is the connection to the random phase approxima- 
tion (RPA). The computation of the one-loop effective action implies that the corresponding 
Green's function is obtained in an RPA-in the diagrammatic picture this corresponds to re- 
summing a chain of one-loop fermion bubbles. This corresponds precisely to the derivation 
of the dressed propagator in the hard thermal loop approximation scheme. Thus, in our 
derivation here, we take the opposite approach to the original HTL derivation; we first 
calculate the effective action and from that result deduce the Green's function. The re- 
semblance between the RPA and hard thermal loop approximation also suggests a natural 
direction for improving the effective theory further, and hopefully the world-line method 
will prove to be useful in this well. 

Another interesting extension would be to include also the effects of a finite density in 
our results. This would first of all generalize the derivation of the HTL effective action. 
Secondly, it would be important to study to what extent the non-linear terms cancel also 
at finite densities, since this situation is more likely to resemble situations of astrophysical 
interest . 

Finally, a future topic is how to include the non-static higher order terms. Even though 
this seems to be a formidable task in general, obtaining just the next order correction would 
be of great interest. In particular, a non-Abelian generalization could provide some insights 
into the expansion of the Wilson line in real time P5| , ^ . Such a continuation is not merely 



of academic interest, but could in fact have important consequences for the observables of 



the heavy ion collisions at RHIC ||45[| 
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APPENDIX A: SOME USEFUL IDENTITIES 

In this appendix we give some of the intermediate stages in the calculation of the hard 
thermal loop effective action. With the Fourier expansion 

A,{x,x,) =TY.j 7|^A,(a;„,p-)e-'(----^^-^) , (A.l) 
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where uom = ^mnT is the Matsubara frequency, and using the change of variables in Eq. 



(p!3|), we find the quadratic term in Eq. (|l^) to be, before any gauge-fixing 



m I •' ) •' l^^J "'O n=-oo 



M4(Ti)+ Z4 4 



expjip ■ [u{ti) + z\ +iq- [u{t2) + z\- iu. 

^^A,{u;^,p)A^i-uj,^,-p) Y: (-1)"/ - 



U4{t2)+ Z4 4- 



X 

■nl3T2 



t 



t fTl 



Jo 



M dn dT2 / Vu M^(ri) + 



exp < -it^-m [M4(ri) - U4{t2)] + ■ [M(ri) - u{r2)] - iuJmnl3- 



(A.2) 



We now use the trivial identity uiji) = J dT[u{T)S{T — Xj)], as well as the generalization of 
Eq. ®, 



7- ; A* 



linear in each A 



(A.3) 



with the notation A^-^^ = A(tOrn,p), A^"^^ = A{—iJm, — p), to write Eq. ( |A.2[ ) as 
^soa(2) = -2e^rE / E (-1)" r ^e-"^/^^/- f dr, T ^r, x 



exp 



(2vr)3 „i- 
n/5(4'^ + Af ) _ .u^n/Sin - T2) 



M / Vu exp 



t t 
where the source term = ( J^, J4) is given by 



u 



linear in 



i. . (A.4) 



each A 



1 = 1 



Mr) = luon [5{t - n) - 5(r - T2)] - E 4^(5(r - Tk)dr, • 

We can then change the variable of integration, 

Uf, = Uf, -2 dfGBir, f) J^(f) , 



(A.5) 
(A.6) 



(A.7) 



where Gb{ti,T2) is the one-dimensional Green's function on a circle of circumference t, 

in - r2? 



Gb{ti, T2) = \ti - T2\ 



(A.8) 
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and 



dr^GBin, T2) = e(ri - T2) 



2(ri - Ts) 



dr^dr^GB{Ti,T2) = 2 



-6{ti- T2) 



(A.9) 



The resulting path integral then becomes Gaussian, and with the normalization in Eq. 
we obtain (with the implicit condition that only terms containing one power of each A^^^ is 
to be kept), 



soft (2) 



[2tcY Jo t-^ Jo Jo 



°° dt 



87r2" 



{2tt 



exp 



^2 



.UmuPin - T2) , n/3(4'^+Af)' 



t 



exp 


[-1 







"2 X 



rt rt 

dr dfJf,{T)GB{T,f)J^{f] 
10 Jo 



871- 



d^p 



exp 



f"°o dt 



72 X 



JO t'^ JO JO 

P^GBfri, rs) + idrMBiri, r2)x 



exp 



{pk{A^k^ + 4'^) - P4(A« + Af)} - A^^AfdrA.GB{ru T2 



(A.IO) 



where we have put = Um and = p^ + p\ in the last expression. By finally using the 
identity 



/ dTi I dT2f{Ti — T2) = [ dx{l — x)f{xt) 
Jo Jo Jo 



(A.ll) 



we can write Eq. (|A.1CI| ) as 

.2 ^ 



soft (2) 



exp 



rp.n^x + .e,4 {j>,(A« + 4^))} (1 - 2x) - ^^^'^^'i^ ' ^(^)} 

(A.12) 



with e^4 = — 1 if /i = 4 and 1 otherwise. The dependence on the vector field can then be 
expanded again, to recover the specific quadratic terms discussed in more detail below. 



1. The ^4-^4 part 

Consider first the terms involving the A4 field. Since we are using the gauge condi- 
tion ^4^4 = to simplify the actual calculation, Eq. (|A.12| ) changes in that there is no 



P4-dependence when the spatial components Aj. = 0. Instead of a sum over Matsubara 
frequencies the integration over Z4 is left over; in other words, we do a three-dimensional 
Fourier transform of A4 instead of Eq. ( |A.1| ). By expanding the y44-dependent exponential. 
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nP(A^l'>+Af^)/t-^Ai'-'> Af\l~5{x)) 



linear in 



V/?^_2(l-5(x))\ (,) 
- /I4 /I4 , 



each A 



^2 



t 



(A.13) 



and using the following integral representation of the Bessel function of the second kind, 



dt 

2 7o F+i 



exp 



X z 

-2 * + T 



(A. 14) 



we obtain Eq. 



2. The Ak-A.-, part 



Turning now to the spatial parts of the vector field, we first use Eq. ( |A.12| ) again and 
put ^4 = 0. We then use the expansion 



iVh (^1^' ) (l-2aO-2A(i) • A(2) {\-&{x)) It 



linear in 



each A 



kj- 



2(1 



+ (1 - 2x)' 



.(1).{2) 



(A.15) 



to obtain Eq. (^11) . Using Eq. (|A. 14|) and performing a partial integration in x, we then get 
Eq. (|23|). To find the leading T-dependence in the effective action, only the sums over the 
square-root terms in Eq. (|2^) need to be kept, and by changing m ^ — (m + 1) in the second 
square- root we have. 



L 



(2) 

soft (2) 





(a 




T>0 




i 



p,p, A\^Uf \l-2x) 
-p^x"^ + — 2k^p^)x + k\ 



+ 2 / da 



1(1). i{2) 



P2x(l -x) + [ki - apiY 



(A.16) 



where k^ = {2m + l)7rT. By integrating over a and performing a partial integration over x 
in the second term inside the square brackets, Eq. ( [A.16|) becomes. 



r(2) 

-^soft (2) 



T>0 



\ / m=~oo L I J 

^0 \ . /-p2a;2 + (p2 _ 2A;4p4)x + kl ) 



(A.17) 



To proceed from the above equation, we use that 



TT 



(02 + ^)2 • 



(A.18) 



Since we are discarding all terms but the leading T dependence, we can in fact neglect the 



term pkPjA). Aj in the numerator PkPjAl''Ay>{l-2xy. This can be shown by using Eq. 
( |A.18| ) and performing the x-integration. We can then write Eq. (|A.17|) as. 
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r(2) 

-^soft (2) 



T>0 

2 



(A.19) 



In the first and second term, we write the sum over m as a contour integral and keep only the 
finite T part. For the first term the contribution to the contour integral comes from the cut 
of the square- root, whereas for the second the contour can be closed and the contribution 
evaluated by the residue theorem. Remarkably, the two contributions cancel each other and 
we are left with the last term in Eq. ( A.19| ). 



APPENDIX B: THE HIGHER ORDER CONTRIBUTIONS 



In going from the expressions for /^(j) in Eqs. ^^e results in Eqs. (^21) - (^11), 

we follow the general method described in |jl2|. To be specific, consider the the sum over 



winding modes in Eq. ( ^6|) coming from Eq. 



(B.l) 



n=l 



Using a different integral representation for the Bessel function. 



dt cosh(z/)f:)e 



-X cosh{t) 



(B.2) 



we have 



OO „OQ oo 

^ {-l)''n^-^K2-N+k {nf3p) = dt COsh[(2 -N + k)t] ^ (^_l)n^N-2^-nPpcosh(t) 

n=l n=l 



N~2 



COSh[(2 - + fc)t] ^ ^n-nfSpcosHt) 

[pCOsh(t)]^-2 [dpi ' 



-1 



\ df3 



dt 



n=l 

cosh[(2 - A^ + k)t] 



[cosh(t)] 



N-2 



g/3pcosh(t) _j_ ^ 



(B.3) 



Now, with the relation 



cosh[(2 -N + k)t] = cosh[(A^ - 2 - k)t] = 2l^-2-^l-i[cosh(i)]l^-2-'=l + |iV - 2 - A;| 



X 



M , 

E- 

m=l 



-11 



m 



\N-2-k\-m-l 
m — 1 



21^-2-^1-2— i[cosh(t)]l^-2-fc|- 



2m 



(B.4) 



where the sum over m runs up to the largest integer M < \N — 2 — k\/2, we have for the 
different values of /c. 
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-1 



dt 



\dl3 



N-i 

E 

3=0 



( 1 M 

' +E 



2[cosh(t)]2 - 1 
[cosh(t)] 



Af-2 



+ 5, 



k{N-2) 



[cosh(t)] 



JV-2 



+ 



,(1) 



^ \ [C0sh(t)]j- ' [C0sh(t)]^+2m 



(B.5) 



Here a^^j-* captures the numerical factors in Eq. (|B.4|) . 
To proceed, we now use the Melhn transform, 



dx f{x)x 



s-l 



to write 



re+ioo 

fi^) = 7r~ dsg{s)x~' , 

+ 1)"' = 7^ dsT{s){l - 2'-^)as)x- 



(B.6) 
(B.7) 

(B.8) 



Using this equation we can perform the t-integration (the integral over powers of the hyper- 
bolic cosine function gives ratios of F-functions |3^), and by taking the {N — 2) derivatives 
with respect to f5 and inserting the the prefactors of T and p, we finally obtain Eq. (^). 

In scalar QED, the situation is almost identical, with a few notable differences. Firstly, 
the overall result has to be multiplied by (—1/2), to correct for the Fermi statistics and the 
number of degrees of freedom. Secondly, the sum over the winding modes in Eq. ( |B.3|) has 
no factor (—1)", so instead of (— l/[e^ + l]) one has a Bose-Einstein distribution, (l/[e^ — 1]). 
The function g{s) in the Mellin transform then changes to F(s)C(s), and the vertical contour 
in the complex s-plane must now be along Re(s) = 1 + e. 
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